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Abstract 

The purpose of this paper is to present a quantum statistical theory of 2- 



> 
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f^ , dimensional vortex gas based on the generalized Hamiltonian dynamics re- 



cently developed. The quantized spectrum is evaluated for a pair of vortex on 
the basis of the semiclassical quantization rule. This is used to evaluate the 
partition function for a dilute vortex gas. A remarkable consequence is that 



'^ ' the partition function and related quantities are given in terms of the gener- 

fj \ alized Riemann zeta function. The topological phase transition is naturally 



understood as the pole structure of the zeta function. 
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Introduction 

The study of vortex defects in two dimension has revealed many interesting properties 
which have had a wide variety of apphcations in condensed matter physics. For example 
the superfluid ^He plays an important role in understanding the properties of the physics of 
vortex defects. Apart from the quantum fluid such as ^He, it has long been known that in 
classical hydrodynamics the equation of motion for the center of vortices forms a Hamiltonian 
system, (e.g. Kirchhoff equation) [0] where the fluid under consideration is assumed to be in 
an effectively two-dimensional incompressible fluid. In our recent papers we have developed 
a novel theory of vortex motion occurring in the two-dimensional quantum condensates, 
namely, the superfluid liquid Helium as well as the Heisenberg spin model |0|. This 
theory shows that the vortex motion is described by the generalized Hamiltonian dynamics, 
which indicates the deviation from the conventional canonical form. Here we have used the 
time dependent Landau Ginzburg (LG) equation (action) to derive the effective action for 
the motion of center of vortices. The LG action is described by the complex order parameter 
field which is naturally introduced through the use of the coherent state path integral for 
the quantum condensates. 

On the other hand, it is well known that the assembly of vortices occurring in the XY- 
model and the similar two-dimensional condensates reveal the topological phase transition, 
known as the Kosterlitz Thouless (KT) transition [Q . This shows the occurrence of dissolu- 
tion of vortex-antivortex pairs at some finite temperature. This implies also the vortex pair 
excitations as the fluctuation out of a state with quasi long range order in two-dimensional 
XY and superfluid models. 

The purpose of this paper is to study the quantum statistical mechanics of 2-dimensional 
vortex gas on the basis of the generalized Hamiltonian dynamics that is developed in the 
previous two papers. Specifically we are concerned with the statistical properties of vortex- 
antivortex pairs. We are first concerned with the quantization of a vortex pair by using 
the Bohr-Sommerfeld quantization scheme. By using the resultant quantized spectra for a 
vortex pair, we next evaluate the statistical partition function for a vortex pair in a dilute 



gas limit. We get a remarkable consequence that the partition function is given in terms of 
the generalized Riemann zeta function. As the characteristic property of the zeta function 
we see that the topological phase transition of KT type is naturally explained by the pole 
structure of the zeta function. 

Basic Formulation 

We start with a concise review of the previous results for the later discussion. The basic 
standpoint is that the center of vortices (xj, i/i) = Xi are regarded as canonical conjugate 
like variables. The effective action for the center of vortices is given by 

Here the first term is called the canonical term and the coefficients Ai, Bi are given in terms 
of the "density function", say cr, together with the function Oj = log |x — Xj|, 0) whereas the 
second term represents the effective Hamiltonian for the motion of vortex centers which has 
the well known log potential form: 

-f^eff = -^PO X! P«Pi 1°S \^i - Xj\+C (2) 

and hence the equation of motion for the vortex centers becomes 

^GijXj = — —J- (3) 

j oXi 

where G is the metric tensor and /x denotes the strength (charge) of the vortex. C represents 

the chemical potential or self-energy that is given by a sum of the self-energy of isolated 

vortices. The tensor G coincides with the symplectic form (two-form) and it is derived from 

the effective action uj 

duj = J2 GijdXi A dXj (4) 



^ The definition of the density function a is different for two cases of bose and spin fluids; see ref. 

[2] and [3]. 



with 

G= (5) 

where each submatrix g is given by the density function etc. through the coefficients Aj, Bi 
in (|l|) and the concrete form for these is given in term of the overlap integrals between the 
gradient of a as well as the a functions as has been done in the previous papers for the bose 
fluid and Heisenberg spin model. Thus the physical meaning of these terms suggests the 
effect of finite size of vortices. We see that 

^11 = ^^ = 0, g''=g' + Ag. (6) 

gO represents the diagonal and Ag the off-diagonal term of the matrix g^"^ respectively with 
respect to the vortex indices. We can introduce the poisson bracket using the tensor G: 

{A5}PB = E(|^||-{^-5})- (7) 



jj dXi dXj 



As a special case we have 



{X„X,} = (G-%. (8) 
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In the more concrete form 

{x„ |/,}PB = {g'T' = {9%' + {{9')-'^9h + ■■■■ (9) 

The second term of LHS can be understood to mean that Ag causes the anomalous effect. 
Note that the above formulation can be regarded as an important generalization of the 
classical Kirchhoff equation x = {x, if }pb and y = —{y,H}pB , with the Hamiltonian 
has the same logarithmic form, whereas {A, -B}pb denotes the usual Poisson bracket. The 
quantization is carried out by simply replacing the generalized PB by the commutator: 
{xi,yj} -^ [xi,yj]. Furthermore to be mentioned is the following fact: The equation of 
motion for the case that Ag vanishes is rewritten as 

{goh (k X X,) = -f. (10) 

oXi 



This equation suggests that the canonical equation of motion (Kirchhoff equation) is alter- 
natively regarded as the balance equation of two type "forces" , the one is the right hand 
side, that represents the force derived from the "potential function" Hes and the other is 
the left hand side, that represents the velocity of the i-th vortex in a literal sense times the 
vorticity vector "(^o^" , where k means the unit vector perpendicular to (x, y) plane. The 
latter corresponds to the so-called Magnus force, f^. 

Quantized Spectra for a Vortex pair 

We shall now treat the system of two vortices whose charges are opposite each other, 
that is, /ii = —fi2 = ^ and quantize it. The most interesting quantity is the quantized 
energy spectrum for the vortex pair, which will be used to evaluate the statistical function 
of the vortex gas. In order to carry out this, we first consider the case that the non-canonical 
effect is omitted and its effect will be taken into account as a perturbation. By having this 
in mind, we shall introduce the "center of mass" and "relative" coordinates as 

xi+X2 = 2Q, yi + y2 = 2P (11) 

and 

Xi-X2 = q, yi-y2=p (12) 

where (xj, yi) is the coordinate of the center of i-th vortex. The effective Hamiltonian which 
is derivatived previously becomes 

H,s = E= i/iV2 log(g' +p') + C , (13) 

where the Hamiltonian is changed to positive definite. Hence we obtain the equations of 
motion as follows: 



^ It should be noted that this effective action do not have the cut-off parameter. 



together with 

dq pofi'^p dp pofi'^q 



(15) 



dt 2 (g2 + p2) ' dt 2 (g2 + p2) ' 

From the first set of equations, it follows that Q and P become constant of the motion, 
which means that the center of vortices is at rest. On the other hand, from the second set 
of equation it follows that the orbit in the (p,q) plane forms a circle: p^ + g^ = constant. 
We shall carry out by using Bohr- Sommerf eld quatization, namely, 

f^uo = {n + ^)h. (16) 

where n = integer and the factor the factor | is added, which corresponds to the zero point 
energy. f\ Physically speaking, the above quantization rule suggests the quantization of the 
angular momentum carried by a vortex pair. In this way, if we note that the unperturbed 
orbit C for the relative motion becomes a circle of radius , say, R, the quantization for 
non-perturbed case is calculated as 

3-^J JdqAdp={n+^-)h (17) 

which leads to 

^rf = (n + l)/^. (18) 

Hence the radius R takes a quantized value; i?^ = -^n with n integers. In the above 
expression, qq takes go = mpofi (bose fiuids), Qq = Jfih (spin fiuid). Next, if the off-diagonal 
term is taken account of as perturbation, 

Rl, = ^inh-TiR^,)). (19) 

Here the value of F can be written as the function of R. Qualitatively, this becomes small as 
the size of the vortex small and it tends to a constant value as ^ oo (the explicit evaluation 



^ The more exotic quantum condition may be obtained, which results in the quantum number of 
quater integer, see 0, but such a possibiUty is not considered hereafter. 
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may be given in the separate paper) 0. In this way, Rn' should be determined from the 
self-consistent equation, in other words, the quantum number n changes to the modified 
value n', which may be written as 

n^n' = n-'^. (20) 

h 

If we assume that the effect of F is sufficiently small and hence put r(i?„/) = — F = constant, 
then we have the simpler form for the quantum condition: 

Rl' = ^iin + hh + f) (21) 



which yields 



r 

n' = n+-. (22) 



If we substitute (^2]) into (0), we get the energy spectra for the relative motion: 

En = -fx'^Po\ogRl, + C 

= l^'po log(^((n + hh + f )) + C. (23) 

Summarizing the above result, the angular momentum quantum number does not take 
integer values. 



Partition Function of Dilute Vortex Gas 

Now we investigate statistical mechanics of two-dimensional gas of vortices with strength 
/i (charge) in the medium of "charge" neutrality and each vortex interacts through log 
potential. The thermodynamical quantities are obtained through the partition function 
Z = Ti e'^^"^ , where (3 = 1/kT and T is the temperature and k denotes the Boltzmann 
factor. The trace means the sum over all possible states the system is allowed to take. 
Here we restrict ourselves to a special situation: namely, the chemical potential C, which is 
necessary to create a pair of vortex of plus and minus charge, is very large. In this extreme 
situation, the vortex system can be treated as a dilute gas consisting vortex pairs. In this 



limit, it suffices to treat only one pair of vortices to consider the partition function. This is 
easily evaluated to be 



Z = £exp[-/5/7, 



eff 



n=0 



oo 



^exp 



n=0 

which leads to 



-/3|^/.Volog(^((n + ^)/. + f)) + C 



(24) 






'-[^) C(^.5 + j)exp[-/3Cl (25) 

and the free energy of the system is given by 

F = -kT log Z. (26) 

Here C('5, a) is defined by 

oo -r 

which is just a generalized form of Riemann zeta function (in mathematical literature, this 
is called "Hurwitz zeta function"). |^ The conventional zeta function is obtained by setting 
r = 0. Remarkable is that ( diverges at the point at s = 1; namely, ( function has a 
simple pole at s = 1. This pole singularity causes the divergence of the free energy at 
finite temperature; ^^^^° = 1, namely, the critical temperature is given by T^ = ^^^ This 
fact suggests that a phase transition occurs at this temperature. We shall also examine the 
susceptibility due to these pair, which is given by 

X = l^i'P <R'> (28) 

where < R^ > represents the mean square separation of the each dipole pair, which becomes 



<i?2^ =^TT{R'exp[-f3H,s]} 

/•/^mVo _ 1 i I In 



^^ 4 ' 2 "•" /i^ 



(29) 



The above expression shows up the following remarkable consequence; the pole of (p9D, 
namely, 4^° — 1 = 1, defines another transition temperature: If the temperature approaches 
to Tc from the below, the mean separation between the pair diverges at this temperature. We 
see that this position corresponds to nothing but the temperature T^ of the KT transition at 
which a dipole pair is made to dissociate. Note that this temperature is just half of T^ noted 
in the above. The free energy or partition function itself remains finite at this temperature, 
namely, this feature is seen from the formula 

Z = C(2.i + ^) (30) 

which is finite. In the above argument we assume that the position of pole and hence the 
transition temperature is not affected by the presence of shift of quantum number by an 
amount of F, which follows from the theorem concerning the Hurwitz zeta function. 

Summary 

To summarize the present result is that after having calculated the energy spectra in- 
cluding the effect of the non-canonical term coming from the specific nature of quantum 
vortices, we have shown that the topological phase transition of KT type is characterized by 
the pole structure of the generalized Riemann zeta function. Up to our knowledge, this fact 
has not been known by anyone previously, which suggests that characteristics of statistical 
thermodynamics of vortex gas is made clear by exploiting the specific properties of zeta func- 
tion. In this connection, we note the appearance of two different transition temperatures; 
Tc and T^. It seems puzzled to have two different transition temperatures, but it should be 
noted that the present theory is based a rather naive assumption, that is, only one vortex 
pair is considered and no correlation, which arises from the other vortex pairs, is not taken 
into account. If we incorporate this correlation in a proper manner, we may have a correct 
value of the temperature for the topological phase transition. As the other problems we 
must consider the case that the noncanonical term has dependence of the quantum number 
n which is shown by the self-consistent nature of the Bohr-Sommerfeld equation. The details 



of these discussion will be given elsewhere. 
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